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We study the tetraquark states with I G J PC = 1~ 1 h in the QCD sum rule. After exhausting all 
possible flavor structures, we analyses both the SVZ and finite energy sum rules. Both approaches 
lead to a mass around 1.6 GeV for the state with the quark contents qqqq, and around 2.0 GeV for 
the state with the quark contents qsqs. The flavor structure (3 ® 6) © (6 ® 3) is preferred. Our 
analysis strongly indicates that both 7ri(1600) and 7ri(2015) are also compatible with the exotic 
tetraquark interpretation, which are sometimes labeled as candidates of the 1 h hybrid mesons. 
Moreover one of their dominant decay modes is a pair of axial- vector and pseudoscalar mesons such 
0© ■ as 6i(1235)7r, which is sometimes considered as the characteristic decay mode of the hybrid mesons. 
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(*C) ' I. INTRODUCTION 

Hadrons beyond the conventional quark model have been studied for more than thirties years. For example, Jaffe 
suggested the low-lying scalar mesons as good candidates of tetraquark states composed of strongly correlated diquarks 
CIh' in 1976 [![. Especially there may exist some low-lying exotic mesons with quantum numbers such as {J PC ) = (1 '") 
which qq mesons can not access @, y|. However the hybrid mesons with explicit glue can carry such quantum numbers. 
• The experimental establishment of these states is a direct proof of the glue degree of freedom in the low energy sector 
C"| \ of QCD and of fundamental importance. 

The mass of the non-strange exotic hybrid meson from lattice QCD simulations includes: 2GeV [J], 1.74 GeV [f|, 
and 1.8 GeV [f3]. The mass of its strange partner is 1.92 GeV Q and 2 GeV [f|. The hybrid meson mass from the 
constituent glue model is 2 GeV 0] while the value from the flux tube model is around 1.9 GeV @,[9[. The prediction 
from the QCD sum rule approach is around 1.6 GeV [T^.[Tl[. However. Yang obtained a surprisingly low mass around 
Q\ , 1.26 GeV for the 1 h hybrid meson using QCD sum rule 

0\ ■ Up to now, there are several candidates of the exotic mesons with I G {J PC ) = 1 (1 + ) experimentally. They are 
^ i 7^(1400), 7Ti(16G0) and tti (2015). Their masses and widths are (1376 ± 17, 300 ± 40) MeV, (1653±Jf , 225^) MeV, 
(2014 ± 20 ± 16, 230 ± 21 ± 73]MeV, respectively [3 . tti(1400) was observed in the reactions n-p -> rjTT°n [II]; 
pp — > 7r°7r ry and pn — > 7r~7r°?7 ir~v — > rjTT~p [la ]. 7Ti(1600) was observed in the reaction ir~p — » r]'ir~p (r)' decays 
to r)ir + n~ with a fraction 44.5%) [17]]. Both 7Ti(1600) and 7Ti(2015) were observed in the reactions ir~p — » ujTt~ir°p [l8j 
and ir~p — > rjir + ir~ ir~ p [19]. However, a more recent analysis of a higher statistics sample from E852 37r data found 
no evidence of 7Ti(1600) [201 ]. All the above observations were from hadron-production experiments. 

Recently, the CLAS Collaboration performed a photo- production experiment to search for the 1 ^ hybrid meson 
in the speculated 37r final state in the charge exchange reaction ■yp — > 7r + 7r + 7r~ (n) [21]. If 7Ti(1600) was an hybrid 
state, it was expected to be produced with a strength near or much larger than 10% of the 02(1320) meson from 
the theoretical models 22]. However 7Ti(1600) was not observed with the expected strength. In fact its production 
rate is less than 2% of the 02(1320) meson. If the 7Ti(1600) signal from the hadron-production experiments is not an 
artifact, the negative result of the photo-production experiment suggests (1) either theoretical production rates are 
overestimated significantly or (2) 7Ti(1600) is a meson with a different inner structure instead of a hybrid state. 

In fact, the tetraquark states can also carry the exotic quantum numbers I G (J PC ) = 1~(1 It is important 
to note that the gluon inside the hybrid meson can easily split into a pair of qq. Therefore tetraquarks can always 
have the same quantum numbers as the hybrid mesons, including the exotic ones. Discovery of hadron candidates 
with J PC = 1 does not ensure that it is an exotic hybrid meson. One has to exclude the other possibilities 
including tetraquarks based on its mass, decay width and decay patterns etc. This argument holds for all these 
claimed candidates of the hybrid meson. 
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Tetraquark states in general have a richer internal structure than ordinary qq states. For instance, a pair of quarks 
can be in channels which can not be allowed in the ordinary hadrons. The richness of the structure introduces 
complication in theoretical studies. Therefore, one usually assumed one or a few particular configurations which are 
motivated by some intuitions. 

Recently, we have developed a systematic method for the study of multiquark states in the QCD sum rule, and 
particular applications have been made for several tetraquark states [111, HH EH • Our method is essentially based on 
complete classification of independent currents. By making suitable linear combinations of the independent currents 
we can perform advanced analysis as compared with the analysis of using only one type of current which limits the 
potential of the OPE, and sometimes leads to unphysical results. 

In this paper, we first classify the flavor structure of four-quark system with quantum numbers J PC = 1 h . We 
find that there are five iso- vector states. Then we construct tetraquark interpolating currents by using both diquark- 
antidiquark construction ((qq)(qq)) and quark-antiquark pairs ((qq)(qq))- We verify that they are just different bases 
and can be related to each other. Therefore they lead to the same results. By using diquark-antidiquark currents, 
we perform the QCD sum rule analysis, and calculate their masses. Our results suggest that 7Ti(1400) may not 
be explained by just using tetraquark structure, and 7ri(1600) and 7ri(2015) could be explained by the tetraquark 
mesons with quark contents (qq)(qq) and (qs)(qs) respectively. The diquark and antidiquark inside have a mixed 
flavor structure (3 ® 6) © (6 <8> 3). 

This paper is organized as follows. In Sec. [Til we construct the tetraquark currents using both diquark (qq) and 
antidiquark (qq) currents. The tetraquark currents constructed by using quark-antiquark (qq) pairs are shown in 
Appendix [XJ In Sec. IIII1 we perform a QCD sum rule analysis by using these currents, and calculate their OPEs. In 
Sec. IIVI the numerical result is obtained for their masses. In Sec. [V] we use finite energy sum rule to calculate their 
masses again. We discuss the decay patterns of these 1 h tetraquark states in Sec. EU Sec. I VIII is a summary. 



In order to construct proper tetraquark currents, let us start with the consideration of the charge-conjugation 
symmetry. The charge-conjugation transformation changes diquarks into antidiquarks, while it maintains their flavor 
structures. If a tetraquark state has a definite charge-conjugation parity, either positive or negative, the internal 
diquark (qq) and antidiquark (qq) must have the same flavor symmetry, which is either symmetric flavor structure 
6f (g)6f (S) or antisymmetric flavor structure 3f ® 3f (A), and can not have mixed flavor symmetry neither 3f ®6f nor 
6f ® 3f (M). However, combinations of 3f © 6f and 6f © 3f can have a definite charge-conjugation parity. Therefore, 
in order to study the tetraquark state of I G J PC = 1~1 h , we need to consider the following structures of currents 



where q represents an up or down quark, and s represents a strange quark. The flavor structures are shown in Fig. [JJ 
in terms of SU (3) weight diagrams. The quark contents indicated at vertices follow the ideal mixing scheme for inner 
vertices where the mixing is allowed. In the SU (3) limit, the quark contents are suitable combinations of the ones 
shown in this figures. However, the strange quark has a significantly larger mass than up and down quarks (current 
quark mass), and so, the ideal mixing is expected to work well for hadrons except for pseudoscalar mesons. The flavor 
structure in the ideal mixing is also simpler than that in the SU(3) limit. Therefore, we will use the ideal mixing in 
our QCD sum rule studies. 

In the following subsections, we first construct currents by using diquark (qq) and antidiquark (qq) currents, and 
then we show the currents with explicit quark contents. The currents constructed by using quark-antiquark (qq) pairs 
can be related to these diquark currents, and are shown in the Appendix. [XI The tensor currents rj^ v (n^ — —T) VIJl ) 
can also have I G J PC = 1~1 By using tensor currents, we obtain the similar results, which will be shown in our 
future work. 



We attempt to construct the tetraquark currents using diquark (qq) and antidiquark (qq) currents. For each state 
having the symmetric flavor structure 6f ©6f (S), there are two (qq)(qq) currents of J PC = 1 h , which are independent 



II. TETRAQUARK CURRENTS 



<7<7<7<?(S) ,qsqs(S) 
qsqs(A) 
qqqq(M) , qsqs(M) 



6 f © 6 f (S) , 
3 f © 3 f (A) , 
(3 f © 6 f ) © (6 f © 3 f ) (M) 



A. (qq)(qq) Currents 




( lTa C l5q2b(q3a1 l il5CqJ b + qzbl )J.lZ C <fia) + Qla C 1^l5q2b(q3aj5CqJ b + q3blbCqJ a ) , 
qTa C l V q2b(q3aV^Cql b - q 3 bO-^CqJ a ) + ql a Co-^q 2b (q3aYCqJ b - q 3 bl"CqJ a ) , 



(1) 
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(qq)®(qq)= 3 f ®3 f (A) 8 



/ 

udus udds 
dsns/ dsds \usds f 

3§J_ X X 

udud 



dsud usud 




ds(uu-dd) us(uu-dd) u ^d ds(uu+dd) us(uu+dd) dsuu 



FIG. 1: Weight diagrams for 6f ® 6f (S) (top panel), 3f ® 3f (A) (middle panel), and 3f ® 6f (M) (bottom panel). The weight 
diagram for 6f ® 3f (M) is the charge-conjugation transformation of the bottom one. 



where the sum over repeated indices (//, v, ■ ■ ■ for Dirac spinor indices, and a, b, ■ ■ ■ for color indices) is taken. C is 
the charge-conjugation matrix, q\ and q 2 represent quarks, and 53 and (74 represent antiquarks. For the antisymmetry 
flavor structure 3f ® 3f (A), we also find that there are two independent (qq)(qq) currents, 

= Qi a Cl5l2b{ciaalnl5Cql b - mi^sCqla) +ql a Cw 5 q 2 b(qa a 75C$ b - q 3b l5CqJ a ) , (2) 
^4 = qIa C l U Q2b(q3a(J^CqJ b + q 3 b<J^CqJ a ) + q[ a Ca tlv q 2b (qs a 'r' / CqJ b + qsbl u CqJ a ) , 

For each state containing diquark and antidiquark having either the flavor structure 3f ® 6 f or 6f ® 3f , there are no 
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currents of quantum numbers J PC = 1 + . However, their combinations (3f ® 6f) © (6f (g> 3f) can have the quantum 
numbers J PC = 1 h . We first define the currents fp^ L which belong to the flavor representation 3f (g> 6f , and the 
currents ipf^ R which belong to the flavor representation 6f ® 3f separately. We find the following four independent 
currents: 



ML 



4 /.i 



= qia C(J ^l5q2b{qiaYlbCql b + qzbl v lr-,CqJ a ) , 

= qia C q2b{q?,al^CqJ b - q~zbluCql a ) , 

= qia C l V l5q2b{qz a a l j.vlhCql b - q^^ulbCqJa) ■ 



K R 
K R 
K R 
K R 

They all have quantum numbers J p 
definite charge-conjugation parity, 



qia C( i?b(q3aiv c qAb + qzbi^Cqla) , 

: qia C l U l5q2b{qza^^l5Cql b + q 3 bCr^l5CqJ a ) , 
: qia C 1^q2b(q3aCqJ b - q 3 bCqJ a ) , 

qia C(T ^i^q2b{q 3 ai"ihCql b - mi v ibCql a ) . 

1~ but no good charge-conjugation parity. However, their mixing can have a 



MR 



(3) 



where the + and — combinations correspond to the charge-conjugation parity positive and negative, respectively. In 
the present work, we only consider the positive one. 



B. Iso- Vector Currents 



For the study of the present exotic tetraquark state, we need to construct iso-vector (/ = 1) currents. There are two 
isospin triplets belonging to the flavor representation 6/ (8 6/, one isospin triplet belonging to the flavor representation 
3f ® 3/, and two isospin triplets belonging to the flavor representation (3f (8) 6 A ffi (6f (8) 3f) (Fig. [1]). For each state, 
there are several independent currents. We list them in the following. 

1. For the two isospin triplets belonging to 6/ ® 6/ (S): 

f Via = ^iniqqqq) ~ ulc^d^ua^^cdl + u^^cS^) + u^Cj a j 5 d b (u a ^ 5 CdJ + u b j 5 Cd^) , 

\vC = ^(qqw) ~ "^CY 'd b (u a a ^CdJ - Uba^CdJ) + u^Ca, J/V d b (u a 'y' / CdJ - u b j"Cd%) , 
f vL = ^ia(qsqs) ~ u T a C^s b {u a ^^ 5 Csl + u^^Csl) + ulC-i^r oSb (u a ^ 5 Csl + u b j 5 Csl) , 
I Via = V^(gsp) ~ u^Cj v s b (u a a^ v Csl - Ubcr^Csl) + u^Ca^s^Ua^CsJ - U b j u Cs^) . 

where 77^ and i]^ are the two independent currents containing only light flavors, and and jy^ are the two 
independent ones containing one ss quark pair. 

2. For the isospin triplet belonging to 3/ ® 3/ (A): 

f Via = ^tAqsqs) ~ u^C^Shiua^^Csl - Uhj^CsJ) + u T a ' Cj^ 5 s b (u a ^ 5 Csl - u^Csl) , 
1 vfu = Wuiqsqs) ~ "Icy s b {u a a^Csl + Ubcr^Csl) + u T a Ca^ v s b {u a ^i v Csl + u b j"Csl) , 

where and i]^ are the two independent currents. 

3. For the two isospin triplets belonging to (3/ <g> 6/) © (6/ ® 3/) (M): 

{V™ = ^{qqqq) ~ u^Cj^dhiuaCdJ + UhCdJ) + ulCd b (u a ^CdJ + Uhj^CdJ) , 
= ^il(qqqq) ~ u r Ca tJ ,^ 5 d b (u a ^^ 5 CdJ + Ub^lsCdfc) + u^CYl^dbiuaa^zCdJ + u&o>7 5 C<iJ) , 
v¥a = ^(????) ~ u T a Cd b {u ala CdJ - Ubl^Cd?) + u^C^d b (u a CdJ - u b CdJ) , 
vfu = *Pfu(qqqq) ~ u l c l l 'l5M u aO u^5CdJ - u b a^ 5 C3^) + ulCcj^sd^UaYlsCdJ - u b YlbCdZ) , 

{vil = ^fl(qsqs) - ulCj^ShiuaCsf + u b Csl) + ulCs b {u alll Csl + Uhj^CsJ) , 
Vffu = ^¥a{qsqs) ~ u^Ca ^^s^UaY l^Csl + u^lzCsl) + u rr Cj"'-f 5 s b {u a a^'-f 5 Csl + u b a^ 5 Csl) , 
f]fa = ^(qsqs) ~ ulCshiua-ff.Csl - u b ^ u Cs^) + ulC^f,,s b (u a Csf - u 6 Cs^) , 
fl&n = ^tliqsqs) ~ u rr Cj"'y 5 s b (u a (T^ 5 Csl - uw^^Csl) + ulCa^^s^UaYlhCsl - u b YlhCsl) , 
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where rj\ 2 3 4 are the four independent currents containing only light flavors, and r/f 2 3 4 ar e the four independent 
ones containing one ss quark pair. 

We use ~ to make clear that the quark contents here are not exactly correct. For instance, in the current rj^ , the 
state usus does not have isospin one. The correct quark contents should be (usus — dsds). However, in the following 
QCD sum rule analysis, we shall not include the mass of up and down quarks and choose the same value for (uu) and 
(dd). Therefore, the QCD sum rule results for with quark contents usus and (usus — dsds) are the same. 



III. SVZ SUM RULE 

For the past decades QCD sum rule has proven to be a very powerful and successful non-perturbative method [261.I27 
In sum rule analyses, we consider two-point correlation functions: 



U^(q 2 ) = i J d 4 xe^ (01^(^)^(0)10), (4) 
where rj^ is an interpolating current for the tetraquark. The Lorentz structure can be simplified to be: 

IV (<z 2 ) = (^ - ^)n (1) (? 2 ) + ^n<«V) . (5) 

q A q A 

We compute Tl(q 2 ) in the operator product expansion (OPE) of QCD up to certain order in the expansion, which 
is then matched with a hadronic parametrization to extract information of hadron properties. At the hadron level, 
we express the correlation function in the form of the dispersion relation with a spectral function: 



1*00 

n«(<z 2 )= / 



~^-ds, (6) 

s — q z — is 

where the integration starts from the mass square of all current quarks. The the spectral density p(s) is defined to be 

p(s) = J2 § ( s - M n)(0\v\n)(n\rf\0) 

n 

= fy S ( s - m y) + higher states . (7) 

For the second equation, as usual, we adopt a parametrization of one pole dominance for the ground state Y and a 
continuum contribution. The sum rule analysis is then performed after the Borel transformation of the two expressions 
of the correlation function, ([4| and (J6]) 

/>oo 

n <««)(Mi) = SmjIT^V) = / e- s / M *p{s)ds. (8) 



Assuming the contribution from the continuum states can be approximated well by the spectral density of OPE above 
a threshold value sq (duality), we arrive at the sum rule equation 



so 



H(Ml) = /£ e- M v' M B = / e- s ' M Bp(s)ds. (9) 



Differentiating Eq. ([9]) with respect to X/M% and dividing it by Eq. (J9]), finally we obtain 

Y II(Afi) f^e-°/ M Bp(s)ds ■ 1 j 

In the following, we study both Eqs. ([9]) and (flQ)) as functions of the parameters such as the Borel mass Mb and the 
threshold value so for various combinations of the tetraquark currents. 

We have performed the OPE calculation up to dimension twelve. Here we only show the results for currents r/^ 
and 775 1 , which have quark contents qqqq and qsqs, respectively. Others are shown in the Appendix. [Bj 



nf (Ml) = P 

Jo 



1 4 (g 2 s GG) 2 , (qq)\ , (qq)( 9s qaGq) 



TS -7— — -S + — — rS 



18432tt 6 18432tt 6 18tt 2 12tt 2 

( {gsqo-Gq) 2 _ 5(g 2 GG)(qq) 2 ^ If 
V 48tt 2 864tt 2 / Ml \ 



- s ' M *ds (11) 



(g s qaGq) 2 5(g 2 GG)(qq) 2 \ , 1 / 32.g 2 (gg) 4 , (g 2 GG) (qq)(g s qcrGq) \ 



81 576vr 2 
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tjM( M -2s _ f° 1 A 17m 2 3 / (g|gg) m B (gg) m fl (as) ^ 2 / (gg) 2 (gg)(as) , * 

5 1 s; ~ i 4m2 18432tt6 7680^ 6 M 18432tt 6 96^ 4 + 48tt 4 J + \ 36tt 2 + 9^ 2 1 ' 



(ss) 2 _ m s (g s qaGq) m s (g s scrGs) m 2 (g 2 GG) \ _ (qq)(g s qaGq) (qq)(g s s(jGs) 
36tt 2 48tt 4 96tt 4 4608tt 6 / S 24tt 2 12tt 2 

(ss)(g s qaGq) (ss)(g s saGs) m s (g 2 GG)(qq) m 2 (gg} 2 m 2 (qq)(ss) m 2 (ss) 2 



S / M *ds 



12tt 2 24tt 2 256tt 4 6tt 2 2tt 2 24tt 2 

(g s qaGq) 2 (g s qaGq)(g s saGs) (g s saGs) 2 5(g 2 GG)(qq)(ss) 2m e (qq) 2 (ss) Am s (qq)(ss) 2 



V 96tt 2 24tt 2 96tt 2 864tt 2 3 

5m s (g 2 GG)(g s qaGq) m 2 (ss)(g s qaGq) m 2 (qq)(g s saGs)\ l / 32.g 2 (gg) 2 (ss) 2 



4608tt 4 4tt 2 6tt 2 J ' M 2 \ 81 

(g 2 GG)(qq)(g s saGs) (g 2 GG)(ss)(g s qaGq) _ 2m s (qq) 2 (g s saGs) _ 5m s (qq)(ss)(g s qaGq) 

1152tt 2 1152tt 2 9 9 

m s (qq)(ss)(g s saGs) i m s (ss) 2 (g s qaGq) ^ m 2 s (g s qcrGq) 2 m 2 (g s qaGq)(g s saGs) 



9 9 24tt 2 24tt 2 

In the above equations, (ss) is the dimension 13 = 3 strange quark condensate; (g 2 GG) is a D — 4 gluon condensate; 
(gscrGs) is D — 5 mixed condensate. There are many terms which give minor contributions, such as (g 3 G 3 ), and 
we omit them. As usual, we assume the vacuum saturation for higher dimensional condensates such as (0|gggg|0) ~ 
(0|gg|0) (0 1 gg 1 0). To obtain these results, we keep the terms of order 0(m 2 ) in the propagators of a massive quark in 
the presence of quark and gluon condensates: 

iS ab = (0|T[g a (x)g b (0)]|0) 

S ab x 2 m q S ab i6 ab m q (qq) . ^ ab m 2 . 

+ W (. 9c ^Gg) - ^ + 48 » + ■ (13) 



IV. NUMERICAL ANALYSIS 



In our numerical analysis, we use the following values for various condensates and m s at 1 GeV and a s at 1.7 GeV 

(gg) = -(0.240 GeV) 3 , 

(ss) = -(0.8 ± 0.1) x (0.240 GeV) 3 , 

(g 2 GG) = (0.48 ±0.14) GeV 4 , 

{gsqaGq) = —Mq x (gg) , (14) 
M 2 = (0.8 ± 0.2) GeV 2 , 
m s (l GeV) = 125 ±20 MeV, 
a s (1.7GeV) = 0.328 ± 0.03 ± 0.025 . 

There is a minus sign in the definition of the mixed condensate (g s qaGq), which is different from that used in some 
other QCD sum rule studies. This difference just comes from the definition of coupling constant g s [28l. |34||. 

For the currents which belong to the flavor representation 6f ® 6f (S) and 3f ® 3f (A), the spectral densities turn 
out to be negative in the energy region 1 GeV ~ 2 GeV as shown in Fig. O The spectral densities of these currents 
become positive in the region s > 4 GeV 2 . They may couple to the state 7Ti(2015). However, after performing the sum 
rule calculation, we find that the mass obtained from the currents r]f^ and rjf is larger than 2.5 GeV, for instance, 
we show the mass calculated from the current 77^ in Fig. SJ The curves are obtained by setting M| = 2 GeV 2 (solid 
line), 3 GeV 2 (short-dashed line) and 4 GeV 2 (long-dashed line). The left curves (disconnected from the right part) 
are obtained from a negative Borel transformed correlation function, and have no physical meaning. Therefore, our 
QCD sum rule analysis does not support 7Ti(1400), 7Ti(1600) and 7i"i(2015) as tetraquark states with a flavor structure 
either 6 f <X> 6 f or 3f ® 3f . 

When using the currents rjf^ , the spectral densities are positive as shown in Fig. [3J And so we shall use these currents 
to perform a QCD sum rule analysis. First we need to study the convergence of the OPE. The Borel transformed 
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s [GeV z ] 



s [GeV 2 ] 



FIG. 3: Spectral densities for the current r\ itL . The spectral densities for the currents with the quark contents qqqq are shown 
in the left hand side, and those with the quark contents qsqs are shown in the right hand side. The labels besides the lines 
indicate the suffix i of the current r/f* (£ = 1, ■ • • , 8). 




FIG. 4: The mass calculated by using the current rji^, as functions of so in units of GeV. The curves are obtained by setting 
M§ = 2 GeV 2 (solid line), 3 GeV 2 (short-dashed line) and 4 GeV 2 (long-dashed line). The left curves (disconnected from the 
right part) are obtained from a negative correlation function, and have no physical meaning. 
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correlation function of the current 77^ is shown in Fig. [5l when we take sq — 4 GeV 2 . Besides the first term, which is 
the continuum piece, the D=6 and D=8 terms give large contributions. The D=6 terms contain (qq) 2 and the D=8 
terms contain {qq)(g c qaGq), which are the important condensates. We find that the convergence is very good in the 
region of 2 GeV 2 < M§ < 5GeV 2 . Therefore, in this region, OPEs are reliable. 

8F 1 




2 3 4 5 

Borel Mass 2 [GeV 2 ] 



FIG. 5: Various contribution to the correlation function for the current ?7 5m as functions of the Borel mass Mb in units of 
GeV 1 at so = 4 GeV 2 . The labels indicate the dimension up to which the OPE terms are included. 



The mass is calculated by using Eq. (fTO]) . and results are obtained as functions of Borel mass Mb and threshold 
value so- In Figs. [HJ [H] and [5J we show the mass calculated from currents 77^, r^, rj^ and , whose quark contents 
are qqqq. Although these four independent currents look much different, we find that they give a similar result. From 
figures at LHS, we find that the dependence on Borel mass is weak. From figures at RHS where the mass is shown as 
functions of So, we find that there is a mass minimum for all curves where the stability is the best, ft is 1.7 GeV, 1.6 
GeV, 1.6 GeV and 1.7 GeV for four independent currents, respectively. We find that sometimes the threshold values 
become smaller than the mass obtained in the mass minimum region. This is due to the negative part of the spectral 
densities. We also met this in the study of V(2175). See Ref [25j for details. 




1.0 L 



3 4 
Borel Mass 2 [GeV 2 



1.0 



1.0 L 



1.0 



so [GeV" 



FIG. 6: The mass of the state qqqq calculated by using the current 7/jT, , as functions of M B (Left) and so (Right) in units of 
GeV. 



In Figs. [TU1 [Til [T2l and [T51 we show the mass calculated from currents r}§L, r/gf, VtL an d VsLi whose quark contents 
are qsqs. The results are similar as previous four currents. But now the mass obtained is about 0.4 GeV larger than 
the previous ones. The minimum occurs at 2.1 GeV, 2.0 GeV, 1.9 GeV and 2.0 GeV, respectively. 

In a short summary, we have performed a QCD sum rule analysis for qqqq and qsqs. The mass obtained is around 
1.6 GeV and 2.0 GeV, respectively. There are four independent currents for each case, which give a similar results. 
Their mixing would lead to a similar result, too. Compared with the experimental data, they can be used to interpret 
the states 7Ti(160Q) and tti(2015) of I G J PC = These analyses are very similar to our previous paper [25j |. 

where we studied the state V(2175) by using vector tetraquark currents which have quantum numbers J PC = 1" 
and quark contents ssss. 
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FIG. 7: The mass of the state qqqq calculated by using the current r)^, as functions of M B (Left) and so (Right) in units of 
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The pole contribution 

S™e- s l M lp{s)ds 



(15) 



is not large enough for all currents due to the high dimension nature of tetraquark currents. Another reason is that 
these currents have a large coupling to the continuum, which is difficult to be removed. Therefore, we arrive at a 
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stable mass, but with a small pole. To make our analysis more reliable, we go on to use the finite energy sum rule. 



11 



3.0 



^ 2.5 
O 



I 2.0 



1.5 





- 




s =4 



3.0 3.0 



> 



3 4 
BorelMass 2 [GeV 2 



a 



2.0 S 2.0 



1.5 1.5 



1 ' \ 
1 \ \ 


- 


\\ ^^B=2 












M B =4 





3.0 



2.5 



-2.0 



1.5 



so [GeV ] 



FIG. 13: The mass of the state qsqs calculated by using the current 77^, as functions of Mg (Left) and so (Right) in units of 
GeV. 



FINITE ENERGY SUM RULE 



In this section, we use the method of finite energy sum rule (FESR). In order to calculate the mass in the FESR, 
we first define the nth moment by using the spectral function p(s) in Eq. ([7]) 



W(n,s )= / p(s)s n ds. 
Jo 



(16) 



This integral is used for the phenomenological side, while the integral along the circular contour of radius sq on the 
q 2 complex plain should be performed for the theoretical side. 
With the assumption of quark-hadron duality, we obtain 



W(n, s ) 

The mass of the ground state can be obtained as 



Ha 



= W(n, s ) 



OPE 



Mf-(n, s ) 



W{n+l,s ) 
W(n,s ) 



(17) 



(18) 



The spectral functions pf 1 (s) can be drawn from the Borel transformed correlation functions shown in section [LTT1 
The d = 12 terms which are proportional to l/(q 2 ) 2 do not contribute to the function W(n, so) of Eq. ([To]) for n = 0, 
or they have a very small contribution for n = 1, when the theoretical side is computed by the integral over the circle 
of radius sq on the complex q 2 plain. 

The mass is shown as a function of the threshold value so in Fig. Qj] where n is chosen to be 1. We find that there 
is a mass minimum. It is around 1.6 GeV for currents r/f 1 , rj^ , rj^ 1 and 774 1 , whose quark contents are qqqq, while 
it is around 2.0 GeV for currents rj^ 1 , 77^, r^f and r/§ , whose quark contents are qsqs. Here w e ag ain find that the 
threshold values become smaller than the mass obtained in the mass minimum region. See Ref (25| for details. In a 
short summary, we arrive at the same results as the previous SVZ QCD sum rule. 



VI. DECAY PATTERNS OF THE 1 + TETRAQUARK STATES 

In this paper, we have verified that (qq)(qq) construction and (qq)(qq) construction are equivalent (see Appendix [A - ]) , 
and from the second one we can obtain some decay information. The four independent (qq)(qq) currents lead to 
the same mass, and therefore, we shall study the decay patterns from all these currents. We can obtain the 5-wave 
decay patterns straightforwardly: 

1. The current naively falls apart to one scalar meson and one vector meson: 

: n-i(1600)-fO+(ff(600) 1 /o(980)»-) + l"(p(770) > a;(782)...) > (19) 
tti(2000) -> 0+ (cr(600), k(800) •••) + !" 0(770), iT*(892) • • • ) • 
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2. The current naively falls apart to one axial-vector meson and one pseudoscalar meson: 
^ : tti(1600) ->l+(ai(1260), &i(1235) ■■■) + 0" (tt(135) •■•) , 



tti(2000) -> l+(ai (1260), Jfi(1270), 



(20) 



•0 _ (7r(135),if(498)---) . 



3. The current £3^ naively falls apart to one vector meson and one axial-vector meson: 



: 7^(1600) ->1~ (p(770), w(782) •••) + 1 + (oi(1260),6i(1235) •••) , 
tti(2000) -» 1" (p(770), if*(892) • • • ) + 1+ (a x (1260), #i(1270) • • • ) • 



(21) 



4. The current £4^ naively falls apart to one axial-vector meson and one vector meson: 



: tti(1600) -» 1+ (ai (1260), 6i(1235) • 



1" (p(770),w(782) 



(22) 



tti(2000) -> 1+ (ai (1260), K x (1270) • • •) + 1~ (p(770), A - * (892) • ■ • ) • 



7Ti(2000) contains one ss pair, so its final states should also contain one ss pair, and its decay patterns are more 
complicated than 7i"i(1600). We see that the decay modes (|2~T1) and (|22|) are kinematically forbidden (or strongly 
suppressed) due to energy conservation. The decay modes (fl9| are difficult to be observed in the experiments due to 
the large decay width of scalar mesons (a and n). Moreover, the scalar mesons below 1 GeV are sometimes interpreted 
as tetraquark states, and if so, these decay modes should be suppressed due to the extra qq pair (24|. Therefore, the 
decay modes (|20[) are preferred. The 7Ti meson first decays to one axial-vector meson and one pseudoscalar meson. 
Then the axial-vector meson decays into two or more pseudoscalar mesons. However, the second step is a P-wave 
decay. Considering the conservation of G parity, the decay mode ai(1260)-7r is forbidden. One possible decay pattern 
is that 7Ti(1600) first decays to &i(1235)7T, and then decays to w7T7t. 

We can also check the P-wave decay patterns besides S'-wave decay patterns. We find that the current £lf leads to 
a decay mode of two P-wave pseudoscalar mesons by naively relating 57^759 and <9 M 7T 



7^(1600) -> 0- (71,77,77' •• + 0- (77,77,77' •••) , 
7T X (2000) -> 0- (71,77,77'- -0 + 0- (71,77,77' •••) . 



(23) 



Considering the conservation of G parity, decay modes tttt and 7777 etc. are forbidden, and possible decay modes are 
7T?7 and 7T?7' etc. Summarizing the decay patterns, there are two possible decay modes: P-wave many body decay, 
such as W7T7T, and P-wave two body decay, such as 7777 and 7777'. This is partly consistent with the experiments which 
observe tti(1600) and 711 (2015) in the decay modes 7777', unrn and titttttt. However, the experiment has not observe 
them in the final state 7777. Certainly it is desired to study these decay patterns to obtain more information on the 
structure of the 7TiS mesons. 
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VII. SUMMARY 

In this paper we have performed the QCD sum rule analysis of the exotic tetraquark states with I G J PC = 1~1 
The tetraquark currents have rich internal structure. There are several independent currents for a given set of quantum 
numbers. We have classified the complete set of independent currents and constructed the currents in the form of 
either (qq)(qq) or (qq)(qq). As expected, they are shown to be equivalent by having the complete set of independent 
currents. Physically, this seems to make it difficult to draw interpretation of the internal structure such as diquark 
(qq) dominated or meson (qq) dominated ones. Using the complete set of the currents, one can perform an optimal 
analysis of the QCD sum rule. 

Somewhat complicated feature arises from the flavor structure. We have tested all possibilities for the isovector 
1=1 states. In the SU(3) limit, there are three cases of, in the diquark (qq)(qq) construction, 6 <g> 6, 3 ® 3 and 
(3 <g> 6) © (6 ® 3). We find that the former two cases can not result in meaningful sum rule since the spectral functions 
become negative. On the other hand, the mixed case (3® 6) © (6® 3) allows positive OPE with which we can perform 
the QCD sum rule analysis. Actual currents have been constructed in the limit of the ideal mixing where the currents 
are classified by the number of the strange quarks. Hence the quark contents are either qqqq or qsqs. 

We have then performed the SVZ and finite energy sum rules. The resulting masses are around 1.6 GeV for qqqq, 
and around 2.0 GeV for qsqs. The four independent currents lead to the same mass and couple to a single state 
as shown above. Hence one of our main conclusions is that the higher energy states 7Ti(1600) and 7Ti(2015) are well 
compatible with the tetraquark picture in the present QCD sum rule analysis. On the other hand, any combination of 
the independent currents does not seem to couple sufficiently to the lower mass state tti(1400), which was, however, 
described as a hybrid state by K. C. Yang in Ref. [ll|. He obtained a low mass around 1.26 GeV by using the 
renormalization-improved QCD sum rules. The 7Ti(1400) state seems somewhat special, as the experiments show the 
similarity between 7Ti(1600) and 7Ti(2015) as well as the difference between 7Ti(1400) and the above two states, which 
we have discussed in the introduction. 

We have also studied their decay patterns and found that these states can be searched for in the decay mode of 
the axial-vector and pseudoscalar meson pair such as 6 1 (1235)7r, which is sometimes considered as the characteristic 
decay mode of the hybrid mesons. The P-wave modes 7177,717/ are also quite important. 

It is also interesting to study the partners of 7TiS. Especially, we can study the one with quark contents udss, which 
is at the top of the flavor representation 10 (see Fig. [IJ. It has a mass around 2.0 GeV, and the decay modes are 
K + (su)K°(sd) (P-wave) and KKK (P-wave) etc. BESIII will start taking data very soon. The search/identification 
of exotic mesons is one of its important physical goals. Hopefully the dedicated experimental programs on the exotic 
mesons at BESIII and JLAB in the coming years will shed light on their existence, and then their internal structure. 
More work on theoretical side is also needed. We will go on to study other tetraquark candidates. 
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APPENDIX A: (qq)(qq) CURRENTS 

In this appendix, we attempt to construct the tetraquark currents using quark- antiquark (qq) pairs. For each state 
containing diquark and antidiquark having the symmetric flavor 6f © 6f, there are four (qq)(qq) currents: 

tlH = fea7p759la)(g4675926) + (g3a75?la) (Gibl iSfhOlb) + (?3a7M75?2a)(?4&753 , l&) + (<?3a75<?2a) (<?4b7M75<?lb) , 

tin = (<i3a7 u qia){(i4b<r t ivq2b) + {qz a (y^qia){qibi l 'q2b) + {q3aYq2a)(q4bcr^qib) + {q^ao ^q2a){qibi 1 ' 'qu) , 

tin = A Qb A c 4fea7M75<7l&)(<Z4c75<72d) + {q^albqibjiqAcl^l^d) + ((?3a7fu75?26) (<Z4c75?ld) + fea75g2&)(94c7M759ld)} , 

tin = ^abKd{(q3ai u qib)(q4 C cr^q2d) + fea^ ^uqib){qici 1 ' 'qzd) + {qsai v q2b){qic<y^ v qid) + {q^n^b){qAci v qu)} ■ 
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Among these currents, only two are independent. We can verify the following relations 

1 

3^ 



Moreover, they are equivalent to the (qq)(qq) currents 

/ S ^ (S 1 ^ /-S 

— — "2~?i^ + 2^ 2 ^ ' 

For each state containing diquark and antidiquark having the antisymmetric flavor 3f ® 3f , there are also four 
(qq)(qq) currents which are non-zero: 

Zip = (93a7M759la) (94675926 ) + fea75<7la) (9467^75 926) - (93a7^ 7592a) (94675 9l&) - (9~3a7592a) (94&7n759l&) , 

£4 = (93a7 Iy 9la)(946Cr M1 >'726) + (q3a^tu>Qla) {(libl" 0.2b) ~ (9~3a7"92a) (q4b^^qib) ~ (^a^^^a) {q4bY qib) , 

Z^n = ^abKd{(qZal^lhqib){q~4clbq2d) + (q~3aj5qib)(q~4c7 f ij5q2d) - {q~3al nl^b) (?4c755ld) - (93a75926) (94c7/[i759ld)} , 

64^ = A a6 A cci {(g3 a 7 1/ qi b )(g4 c CT MIy (72<i) + (q3a(J^qu)(q4cl ly q2d) - (q~3aYq2b)(q~4cVniyqid) - (q~3aV^q2b)(q4cj"qid)} , 

where two are independent 

£.tp — ~ — g £2^ • 



They are equivalent to the (qq)(qq) currents 



1 A _ _ 1 {A 1 * (A 

H/i — 2^!^ 2 M ' 

/ A ^ /-A ^ >A 

V% — ~"2~~VL M + 2^2^ ' 



For the currents which have a mixed flavor symmetry, we just show the (qq)(qq) currents which belong to the flavor 
representation 3f (g> 6f . 

ill L = (q~3aqia){q4bJnq2b) ~ {q~3al l >qia){q4bq2b) - {q3aq2a)(q4b-f t ,qib) + {q~3alnq2a){q~4bqib) , 

= {q3al^l5qia){q4blbq2b) - (g3a759la)(g467 AI 75926) - (93a7 M 7592a ) (9467591&) + (93a7592a) (9467^75916) , 
£3^ = (93a7"9la) (946^926) - (93a0>9la) (9467^926) ~ (93a7 ly 92a)(946^„9i6) + (q~3a 0>92a) (9467^916) , 

tiff = (93a7 ,y 759la)(946CT My 75926) - (93aO>759la ) (9467^75926) - (93a7"7592a)(94fcO>759lf>) + (93a0> 7592a) (9467^75916) • 

There are also four currents which have a color 8 C (§3 8 C structure, and they can be written as a combination of these 

^andCr 



color l c <g> l c currents. The relations between </>*f L and £^f L are: 



j.ML _ ± C ML , L (ML _,_ ' (ML 1 C ML 

WlH - 4^1^ + 4^2^ + ^?3 M 4?4 M . 

/ML _ 3i M£ 3i ML l f Mi , ^ (ML 

V% — ^ ?lp "+" ^ ?2^ ~+" 4?3ai "t" 4?4^i i 



./ ML _ 1 C ML , - 1 - t Mi , ' (ML , 1 (ML 

— 4?1m ~T 4?2^i 4?3m + 4^4^i ' 

/ML _ 3i M£ 3i ML 1 M£ 1 ML 

^4 M - -J?^ + X^ 2 ^ + I^m I?4p 

and the currents with J PC = 1~+ arc 



We can obtain similar results for £™ R , which belong to the flavor representation 6f <8> 3f can be obtained similarly, 



% = tit L + tit R - (Al) 
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APPENDIX B: TWO-POINT CORRELATION FUNCTIONS 



In this appendix we show the results for the Borel transformed correlation functions as defined in Eq. 
for the currents r/f, r]^ 1 , rj^f , , rjg 1 , v^f and rjgf are indicated by the same upper and lower indices. 
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